Breakdown of the coherent state path integral: 
Two simple examples 
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We show how the time-continuous coherent state path integral breaks down for both the single-site 
Bose-Hubbard model and the spin path integral. Specifically, when the Hamiltonian is quadratic 
in a generator of the algebra used to construct coherent states, the path integral fails to produce 
correct results following from an operator approach. As suggested by previous authors, we note that 
the problems do not arise in the time-discretized version of the path integral. 
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Path integrals are widely known for being an alter- 
nate formulation of quantum mechanics, and appear in 
many textbooks as a useful calculational tool for various 
quantum and statistical mechanical problems (e.g., [3- 
3] ) . From their inception, there has been the problem of 
writing down a path integral for any system that can be 
described by a Hilbert space equipped with a Hamilto- 
nian. One way to approach this problem is with what 
is known now as the generalized coherent state path in- 
tegral [1, [B| which generalizes the coherent state path 
integral for a harmonic oscillator. The key observation 
with path integration [2| is that, given a Hamiltonian 
H, the propagator, e~ ttH , at some time t can be broken 
up into N slices, {e~ ltH / N ) N , and in between each multi- 
plicative term one inserts an (over-)complete set of states 
parametrized by a continuous parameter; generalized co- 
herent states meet that criteria. If we take N — » oo 
we get the time-continuous formulation. Coherent state 
path integrals have become widely and routinely used in 
many areas of physics (see the many papers collected in 
1). 

Glauber coherent states Q are usually understood as 
the most classical states associated with the harmonic 
oscillator. They obey the classical equations of motion 
for a harmonic oscillator and are minimal uncertainty 
states. Perelomov and Gilmore extended the def- 

inition of coherent states to Lie algebras other than the 
Heisenberg algebra (i.e., the harmonic oscillator algebra). 
Since then, these "generalized" coherent states have been 
used in a number of applications (see 1^, U| for more on 
them). In particular, the coherent states form an over- 
complete basis (with a continuous label) and therefore 
admit a resolution of the identity, which is a necessary 
ingredient for the construction of a path integral. For the 
harmonic oscillator, coherent states are represented by a 
complex number a, but for coherent states constructed 
with su(2) (spin), they are points on the Bloch sphere, 
S 2 . 

For the case of the harmonic oscillator, it is commonly 
known that one can easily go between the normal-ordered 
Hamiltonian (all annihilation operators commuted to the 
right) and the coherent state path integral [l| ; this is due 



to the fact that coherent states are eigenvectors of the an- 
nihilation operator. For the general coherent state path 
integral, the "classical" Hamiltonian in the path integral 
is just the expectation value of the quantum Hamilto- 
nian with a coherent state. This prescription results in 
some notable exactly solvable cases, but all such cases in- 
volve non-interacting terms which are essentially linear 
in the algebra generators used to construct the coherent- 
states. When the Hamiltonian involves terms that are 
non-linear in generators (interactions), this prescription 
seems to fail, as this letter demonstrates. 

Path integrals are widely used for developing perturba- 
tive expansions in terms of Feynman diagrams, for non- 
perturbative techniques (e.g., the instanton method), and 
for deriving effective theories @, HJ , but despite the many 
successes of path integrals, they have been on very shaky 
mathematical grounds (for a small "slice" of this history, 
see |12|). In particular, the spin coherent state path in- 
tegral has sometimes produced (quantitatively) incorrect 
results [13-17] unless the time-discretized version is em- 
ployed [lfftUJj' These problems with the time-continuous 
path integral were mostly solved by Stone et al. [li| by 
identifying an anomaly in the fluctuation determinant 
which added an extra phase to the semi-classical propa- 
gator. Kochetov had also found this phase in a general 
context [13] • Furthermore, Pletyukhov [2l| related the 
extra phase in the spin path integral back to Weyl order- 
ing the Hamiltonian in the case of the harmonic oscillator 
(in the simplest case, Weyl ordering corresponds to sym- 
metrically ordering annihilation and creation operators) . 
Additionally, Weyl ordering has been considered in the 
Bose-Hubbard case in [22j]. Unfortunately, this solution 
does not explain the present breakdown under consider- 
ation. 

In this letter, we outline another problem with the 
time-continuous coherent state path integral. This prob- 
lem manifests itself in two simple examples: (i) the spin- 
coherent state path integral and (ii) the harmonic os- 
cillator coherent state path integral (in particular, the 
single-site Bose-Hubbard model). The single-site Bose- 
Hubbard Hamiltonian is a minimal model that demon- 
strates the problem with the normal-ordered path in- 
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tegral. However, the problem itself is more general 
than the toy model considered here and clearly per- 
sists in more complicated models, including lattice Bosc- 
Hubbard models. We use an exact method of calculating 
the partition function, first used for spin (with H = S z ) 
by Cabra et al. 23], and demonstrate that the exact re- 
sult differs from the correct partition function in the cases 
of both normal-ordering of operators (as prescribed by 
most textbooks) and when using Weyl ordering (i.e., it 
cannot be accounted for with the phase anomaly found 
by Solari and Kochetov 18, [2(| and elaborated on by 
Stone et al. [H). 

We begin with the coherent state path integral for spin 
with the standard SU(2) algebra defined on the operators 
{S x , S y , S z } with [SijSj] = iEijkSk, and we define our 
Hilbert space by taking the matrix representation of the 
SU(2) group in (2s + l)-by-(2s + 1) matrices (s being the 
spin of the system). Irrespective of the algebra, we can in 
general define a Hermitian matrix H that acts on states 
in our Hilbert space, and this will be our Hamiltonian. 
Usually, H is polynomial of algebra generators. 

If \s) is the maximal state of S z in our spin-s sys- 
tem, then we can define spin-coherent states as |n) = 
e~ l ^ s 'e~ l9Sy \s) where (0,(f>) are coordinates on the 
sphere S* 2 along the unit vector n (i.e., a point on 
the standard Bloch sphere). These coherent states are 
overcomplete such that 2|±1 J dn |n) (n| = 1 where 
dn = dej) d(cos 9) is the standard measure on S* 2 . Us- 
ing this continuous, overcomplete basis, one can derive 
the standard path integral for the partition function for 
spin from Z = tr e^ l3H in the standard way [l| discussed 
in the introduction: 

Z' = J X>n(r) exp j-J^dr[- (n(r)|d T n(r)) 

+ <n(r)|H|n(r))]}. (1) 

We call the partition function as given by the time- 
continuous path integral Z' in order to distinguish it from 
Z = tr e~@ H since we will find that in general they may 
not agree. The path integral is over all closed paths (since 
it is the parititon function). The first term in the action 
for Eq. (fTJ), (n|9 T n), is the Berry phase term and in (6, </>) 
coordinates — (n|<9 T n) = — si(l — cos9)d T (j). 

In order to highlight what is wrong with a naive use of 
Eq. (fTJ), we employ the method of Cabra et al. 23|. We 
assume (n(r)|JJ|n(r)) = H(cos9(t)) for some function 
H(x) (this is true if and only if H is diagonal). This 
puts the (f> dependence of the action solely in the Berry 
phase term of the action. We then integrate the Berry 
phase term by parts; the boundary term is just A</>(1 — 
cos 9(0)) with A</> = 4>(f3) — (j>(0) ~ 2irk for any integer 
k and cos#(/3) = cos 9(0) since our paths are closed. We 
must sum over the different topological sectors defined 
by the integer k (i.e., how many times <p wraps around 
the sphere). Thus, our only 4> dependence is multiplying 



dc °* e from integrating by parts, and we use standard 

identity for functional integrals / X>0 ^ dr 0(t)/(t) = 
8(f), to get that cos9 must be constant (i.e., dc ° s6> = 0). 



This 5-function allows us to do the path integral over 
2?(cos#), except for the initial value which we call x := 
cos #(0). Taking all of this into account, the path integral 
can then be written as 



Z' = \^ / die 2 "' 8 ' 1 " 1 '"^' 1 '. 

k — — OO ^ 



(2) 



The sum over k can be evaluated as a sum of delta func- 
tions of the form S(s(l — x) — n) for all integers n. Since 
x is in the interval —1 to +1, only finitely many n con- 
tribute (n = to n = 2s to be exact). We can rewrite 
the sum over n as a sum over m := s — n and we get the 
answer (dropping overall constants) 



,-0H{m/s) 



(3) 



Eq. looks very promising, but H(m/s) is not the same 
as (m\H\m). First let us see where it does work. Take 
the simple Hamiltonian H — S z , then (n|i?|n) = scosfl, 
and thus H(x) — sx. This immediately yields 



z 'h=s z - E 



-/3m 



(4) 



and it is easily calculated (in operator language) that 
Z' H=S = Zh=s z ■ The two methods agree for the partic- 
ular Hamiltonian H = S z (the case considered by Cabra 
et al. [23]). On the other hand, if we take H = S z and 
s = 1, we can evaluate (n|5 2 |n) = i (cos 2 9 + l); from 
which we have 



H(x) 



^ 2 + l) 



(5) 



Thus, Z' H=S 



2e 13 + e ^/ 2 , but this conflicts with 

Zh=s 2 — 2e~P + 1 by more than just a multiplicative 
constant. Thus, we have Z' H=S2 ^ Z H= g2 for s = 1, and 
in fact Z' H=S 2 ^ Z H== gi for all s > 1/2. 

Importantly, the two methods agree for any Hamil- 
tonian when s — 1/2. This comes from the fact that 
any (diagonalized) Hamiltonian for a two state sytem 
(s = 1/2) can be written as H — a + bS z (in fact 
H = S z = 1/4), and the above method gives Z' = Z 
when H = a + bS z . 

Also, if we let our Hamiltonian be H — S z /s 2 , then in 
the limit of s 3> 1 Eq. ([3]) reproduces the correct result. 
This is a general result for Hamiltonians that are finite 
polynomials of S z /s, and suggests that "semiclassically" 
(i.e. s tends to infinity), we will still arrive at sensible 
results. 

Agreement can also be forced by considering H(x) = 
x 2 instead of Eq. (O, but this corresponds to replacing S z 
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with (S z ) in the Hamiltoinian instead of just considering 
(H). In the H = S 2 case, it is the difference between con- 
sidering (S 2 ) and (S z ) ; the latter gives correct results. 
There is no a priori reason for this construction. 

To motivate looking for this same issue in a system 
with the Weyl-Heisenberg algebra (i.e., the harmonic 
oscillator algebra), it is known [24[ that one can con- 
tract u(2) (since we constructed our coherent states for 
spins with su(2)) into the Weyl-Heisenberg algebra by 
considering u(2) = span{So, S x , S y , S z } = u(l) ® su(2), 
where we define [S'o,^] = 0. Then define the opera- 
tors J := So, Ji,2(e) := eS y , x , and J 3 (e) := S + e~ 2 S z 
to get the commutation relations [</3,Ji,2] = T^.i, 
[Ji, J2] = —ie 2 J3 + iJo, and [Jo, Ji] = 0. If we let 
e —> 0, we recover exactly the Weyl-Heisenberg algebra: 
f)4 = spanjl, x, p, a^a} with [x,p] = i, [a^a,x] = —ip, 
[a'ajp] = ix. Observe that S 2 is related to a' a in this 
contraction, so we might suspect that terms quadratic in 
oJa might give problems like those found with the spin- 
coherent state path integral. 

A Hamiltonian that uses the Weyl-Heisenberg alge- 
bra to construct its coherent states is the Bose-Hubbard 
model. For a single site, we can write 



H 



U ( 
— n(n 

2 v 



1) 



(G) 



where n = a^a and the a (a^) is the annihilation (cre- 
ation) operator for the algebra [a,a^] = 1. The form 
n(n — 1) — aWaa comes from the normal ordering re- 
quired from a path integral of the form 



we have Z' ~ 1 + e^ 17 / 2 + • ■ ■ , but Z ~ 1 + e» + e 2 ^ 17 + 
• • • . With different asymptotics, Z and Z' are different 
expressions. Note that if we let fj, — > % in Z' , that we 
will get same result. This substitution for fi corresponds 
to replacing n in Eq. ([5]) by (n) = \z\ 2 when writing down 
our action (so instead of (n 2 ), one gets (n) ). 

We now compare this to the semiclassical result. Still 
considering Eq. ©, let change our algebra slightly to 
incorporate a small parameter (akin to the standard h — > 
for normal semiclassics): the representation index; 
see 7 defined in j20|. We note here that different h's 
change the coherent states \z) in the following way: if 
z = -4^>{% + iy), then x — q/c, y — p/d, and h = h/ (cd) 

(and [a, a'] = h). We have used q and p as the standard 
position and momentum for the harmonic oscillator. Up 
until now we have been considering h = 1. 

We can write the propagator between two coherent 
states \zi) and \zf) using a Hubbard-Stratonovich trans- 
formation and the propagator for the harmonic oscillator: 



K{z* f ,z i ;t) = {z f \e- iHT ' h \z i ) (10) 
J diuei^ + ^ T+ i uhT , (11) 



iT 
2nUh 

where we have defined 



* u = z} Zi e^ T + ^ i(N 2 + \z f \ 2 ). (12) 

In path integral notation we can write out the propaga- 
tor [H 







j V 2 z exp < 









— (z*z — z* z) 
2 V ' 



(7) 



We can solve this path integral with the same method 
used to obtain Eq. Q in the spin-coherent state path 
integral. Let z = ^fne %6 , so that the measure becomes 
T> 2 z = VnDQ and the action becomes S = j dr(in9 — 
[in + y« 2 )- Integrating by parts on the nQ term then 
integrating over T>9 will fix n to be constant, and the 
boundary term will fix n to be an integer. Since n is 
radial, it can only be positive so we directly obtain 



Z' = e Mn 



(8) 



71=0 



But this differs from the partition function that we can 
easily calculate in operator language: 



2 — e )»^-fii(n-l 



)/3 



(9) 



77=0 



We see that a similar problem to that of the spin coherent 
state path integral here. To see it explicitly, for U 3> 1, 



K(z},z t ;T)= / V 2 z exp{$[z,z*]//i}, (13) 

Jz(0)=z I 



where $ = T + S, 
1 



r = -[z;.z(T) + z*(o)z / -|z / | 2 -|z / | 2 ], 

S=-[ dt(zz* - z*z)-i [ dt (z\H\z 
2 Jo Jo 



(14) 
(15) 



Performing the standard semiclassical analysis and al- 
gebra (see and [Tj|) the semiclassical propagator 
takes the form 



K sc (zf, Zi, T) 



E 



-) 

hU 



1/2 



ia 2 $. 



h dui 2 



-i/2 



x exp 



+ I(w + fi)T - zA 



(16) 



where the sum is over solutions to the consistency equa- 
tion given by = or u = -Uz* f z i e' l( -' Aj+ ^ T , and 
we have defined A = + 2u))T. The term A comes 
from the fixing of the fluctuation determinant anomaly 
described in detail by Stone et al. [19] for the SU(2) case. 
However, if we try to get Eq. (|16|) by using the method 
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of steepest descent on Eq. (fTTj) with h —> 0, we won't 
get the same result. This is because of what has been 
shown by others 2(| 21 1: that the semiclassics will give 
results consistent with the Weyl ordering of the Hamil- 
tonian (naively ordering all a and a^'s symmetrically). 
The usual normal ordered Hamiltonian takes the form 
(inserting h's) H = —fin + jn(n — h) while the Weyl 
ordered Hamiltonian takes the form (up to a constant) 
Hw — —fin + ¥n(n + h). If we derive Eq. (TTTI) for H\y, 
we will find the the steepest descent exactly agrees with 
Eq. ([To]) just as expected [20I |21|. 

While the semiclassical result is not a new one, it 
shows that the path integral is not dealing with the 
same Hamiltonian. Unfortunately, our exact calculation 
suggests that the path integral is dealing with H' = 
—fin + ^n 2 while semiclassics suggests it is dealing with 
Hy/ = —fin + ^n(n + 1). These two methods differ but 
both are not the Hamiltonian under consideration. In 
the case of the Weyl ordered Hamiltonian, we can write 
our original Hamiltonian in Eq. (|6]) as H = Hw — Un 
which is Weyl ordered (up to a constant). This ordering 
can be used to modify the path integral by an extra term 
—Un. This correction to the path integral suggested by 
Weyl ordering does not fix the exact calculation as can 
be easily shown, but it does motivate an ad hoc correc- 
tion to the path integral to "fix" our exact calculation. 
We use the following action (going back to h = 1): 



S = 



J dt f- M |2| 2 + ||*| 



(17) 



This action is constructed by just changing the operator 
n to a function |z| 2 ; while this gives correct results with 
the method which gives Eq. ([8]). there is no a priori rea- 
son to suspect this of being the action. Similarly, if in 
the spin-coherent state path integral, we replace the op- 
erator S z with its expectation value (n|5 z |n) everywhere, 
we will get the correct result. This means, in particular, 
for H = S\ that instead of (S 2 ) in the spin-path integral 
we have (S z ) 2 . In general, if one substitutes the genera- 
tors of the coherent states in the Hamiltonian with their 
expectation value, one obtains the correct result for Z 
with the methods used to derive Eq. {3} and Eq. ([8]). 

Corrections aside, a simple way to see what has gone 
wrong is to return to Eq. ([5]). This H(x) function can 
not achieve the value 0, but H = S 2 clearly has such an 
eigenvalue. This is due to the fact that for higher dimen- 
sional representations of SU(2) not every eigenvector of 
S z can be rotated into another with a standard SU(2) 
rotation. On the other hand, the coherent states we used 
are a complete set for even higher dimensional represen- 
tations, so in principle, we should not lose any informa- 
tion about the m — state. Continuity in n seems to 
be the culprit: H(x) came from a time discretized form 



(between time slices j and j + 1) (rij + i\S 2 \nj) , and we 
have (nlS^I — n) =0, so (nj + i\S 2 \iij) can attain zero, 
but not for any paths that are "close" to each other (i.e. 
rij « nj_|_i) as the continuous time path integral assumes. 
As such, the discrete time path integral (before a continu- 
ity assumption is imposed) can unambiguously give the 
correct results to a calculation. 

To conclude, in the time-continuous formulation of 
the path integral, neither the action suggested by Weyl- 
ordering nor the action constructed by normal ordering 
give the correct result when evaluating Z via path inte- 
grals. 
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